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ABSTRACT 

The energy spectrum of the cosmic microwave background (CMB) allows constraining 
episodes of energy release in the early Universe. In this paper we revisit and refine the com- 
putations of the cosmological thermalization problem. For this purpose a new code, called 
CosmoTherm, was developed that allows solving the coupled photon-electron Boltzmann 
equation in the expanding, isotropic Universe for small spectral distortion in the CMB. We 
explicitly compute the shape of the spectral distortions caused by energy release due to (i) 
annihilating dark matter; (ii) decaying relict particles; (iii) dissipation of acoustic waves; and 
(iv) quasi-instantaneous heating. We also demonstrate that (v) the continuous interaction of 
CMB photons with adiabatically cooling non-relativistic electrons and baryons causes a neg- 
ative A'-type CMB spectral distortion of Aly/Iy ~ 10 in the GHz spectral band. We solve 
the thermalization problem including improved approximations for the double Compton and 
Bremsstrahlung emissivities, as well as the latest treatment of the cosmological recombina- 
tion process. At redshifts z < 10^ the matter starts to cool significantly below the temperature 
of the CMB so that at very low frequencies free-free absorption alters the shape of primordial 
distortions significantly. In addition, the cooling electrons down-scatter CMB photons intro- 
ducing a small late negative y-type distortion at high frequencies. We also discuss our results 
in the light of the recently proposed CMB experiment Pixie, for which CosmoTherm should 
allow detailed forecasting. Our current computations show that for energy injection because 
of (ii) and (iv) Pixie should allow to improve existing limits, while the CMB distortions caused 
by the other processes seem to remain unobservable with the currently proposed sensitivities 
and spectral bands of Pixie. 
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1 INTRODUCTION 

Investigations of the cosmic microwave background (CMB) tem- 
perature and polarization anisotropies without doubt have al- 
lowed modern cosmology to mature from an order of magnitude 
branch of physics to a precise scientific discipline, were theoreti- 
cal predictions today are challenged by a vast a mount of observa- 
tional evidence. Since its discovery in the 60's l lPenzias & WilsoiJ 
1 1 9651). the full C MB sky has been mapped by Cobe/Dmr 



unprecedented precision, while both ACT (e.g., see iHaiian et al 



2010; Dunklev et al. 2010 : Das et"ai]|201ll) and SPT jLueker et al 



dSmoot et alj Il992h and Wmap ( iBennett et atl l2003l : IPage et al.l 
l2006h . while at small angular scales man y balloon-borne an d 
ground-based CMB experiments like Maxima teanany et al.'2000l). 
Boomeran g fNetterfie ld et al. 2002), Dasi (Halverson et al. 2002), 
Archeops jSenoit et alj|2003l) . Cbi jPearson et al.ll2003l) and Vsa 
jRubino-Martin et al.l l2003l) ~provided important additional lever- 
age, significantly helping to tighten the joint constraints on cosmo- 
logical parameters. Presently PLANCiflis producing CMB data with 
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l20ia : ,Vandeiiinde et al..20IO,) are pushing the frontier of TT CMB 

S lower spectra at small angular scales. In the near futu re SPT poiH 
McMahonetZll2009h and ACTPoi| jNiemack eta"i]|2010l) will 
provide additional small scale ii-mode polarization data, comple- 
menting the polarization power spectra obtained with Planck and 
further increasing the significance of the TT power spectra. Also 
QuiE-^ dOUIET Collaboration et al1l2010h already now is observ- 
ing the CMB in Q and U polarization with the aim to tighten the 
constraints on the potential level of S-modes at multipoles / < 500, 
while SpiDEpjfl and Polarbeae(^ are getting ready to deliver precision 
CMB polarization data in the near future. 

With the new datasets, cosmologists will be able to determine 
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the key cosmological parameters with extremely high precision, 
also making it possible to distinguish between various models of 
inflation by measuring the precise value of the spectral index of 
scalar perturbations, n^, and constraining its possible running, ;3,^,„. 
Furthermore, with the Planck satellite primordial B-mode polariza- 
tion may be di scovered, if the te nsor to scalar ratio, r, is larger than 
a few percent jThe Planck Collab oration 2006.) . This would be the 
smoking gun of inflation and indicate the existence of gravitational 
waves, one of the main scientific goals of Planck, Quiet, Spider, 
PoLARBEAR and other CMB experiments. 

However, the anisotropies are not the only piece of infor- 
mation about the early Universe the CMB offers us. In addition, 
the energy spectrum of the CMB tells the tale about the ther- 
mal history of the Universe at very early times, well before any 
structures had formed and when the baryonic matter and radia- 
tion were tightly coupled via Compton scattering. It is well known 
that any perturbation of the full thermodynamic equilibrium be- 
tween photons and baryons in the ea rly Universe inevitably leads 
to spectral distortions in the C MB /Zeldovich Sl Sunvaev 'l 969l. 

iSunyaev & Zeldov ichI 1970dl3: IZeld ovich et al. 1972; Sun vaevI 
1974 |m arionov & SunvaevI 1 19741) . This for example could be 
cause d by some relea se of energy by decaying relic parti- 
cles (lHu&Silklll993bh or evaporating primordial black holes 
jCarr et al.ll2010l) . and depending on their lifetim^ and the total 
amount of energy that was transferred to the medium, the residual 
CMB spectral distortion will be smaller or larger. 

The process that erases possible spectral distortions of 
the CMB in the early Universe and thereby attempts to 
restores full thermodynamic equilibrium is called thermal- 
ization. In the past it has been studied extensively, both 
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energy release (z > 2 x 10 ), photon production at low frequen- 
cies by double Compton emission and the up-scattering of pho- 
tons to high frequencies by Compton scattering are so efficient that 
the thermalization process is practically perfect for nearly arbitrary 
amounts of ene rgy release, such that no distortion should remain 
today (e.g., see iBurigana et al]|l991bl) . However, below this red- 
shift the CMB spectrum becomes vulnerable, and spectral distor- 
tions that could still be observable today can form. 

In connection with early energy release two types of CMB dis- 
tortions, called II- and y-type distortion s, are usually distinguished 
(e.g., see llllarionov & Sunvaev|[l975al lbl). The first type of distor- 
tion is characterized by a frequency-dependent chemical potential, 
^(v), that is very close to constant at high frequencies, and vanishes 
at very low frequencies. This type of distortion forms at times when 
the Compton process is still able to achieve full kinetic equilibrium 
with the electrons (redshifts z ^ few x 10^), and the CMB spec- 
trum reaches a Bose-Einstein distribution with occupation number 
= l/[e*"+''<'') - 1], where x = hv/kT^, with denoting the CMB 
temperature. However, photon production by double Compton and 
Bremsstrahlung emission stopped being very fast, so that only at 
low frequencies full equilibrium can be restored. 

At low redshifts (z < few x 10' - 10"*), on the other hand, 
y-type disto rtions form. These are also k nown in connection with 
SZ clusters l lZeldovich & Sunvaev|[T969l) and are characterized by 



For black holes the evaporation timescale is t^y oc M jHawkinJl974l) . 



a constant temperature decrement at low and an increment at high 
frequencies. This type of distortion is formed when the Compton 
process ceases to be extremely efficient, so that full kinetic equi- 
librium between electrons and photons can no longer be achieved. 
Photons produced at low frequencies by Bremsstrahlung are there- 
fore not up-scattered very much and remain close to the frequency 
band they were emitted. 

Both t ypes of distortions were tightly co nstrained by 
Cobe/Firas dMather et 'Zl l 19941 ; iFbTsen et all 1 19961) only leaving 
room for |y| < 10"' and |p| <, 10^* at cm and dm wavelength. 
This imp oses strong limits on possible energy release by decaying 
particles i IHu & Silkl 19"93bh . which for some particle lifetimes even 
today are compa rable to those derived, for example, from big ban g 
nucleosynthesis l lKusakabe etalibOOd ; iKohri & Takahashill2010t) . 
However, since the 90's detector technology has advanced very 
much. Already about ten years ago, a Firas type-II experiment 
with roughly 100 times better sensiti vity in principle would have 
been possible jFixsen & Matheill2002l) . Additional new technology 
was used in several flights of the Arcade balloon-borne instrument 
jKogut et al. 2004; Fixsen et al. 2011; Kogut et al. 2011), from 
which improv ed limits on the CM B temperature at v ~ 10 GHz 
were derived dSeiffertetallbOIlh . However, until now no signa- 
tures of primordial distortions were found, and the limits on overall 
CMB spectral distortions derived from Cobe/Firas until today re- 
main most stringent. 

However, very rec ently a new space-based mission, called 
Pixie jKogut etal .11201 ih , was proposed. Pixie will have exquisite 
spectral capabilities with 400 channels at frequencies ranging from 
V = 30 GHz to 6THz. The claimed sensitivities should allow to 
tighten the constraints on yu- and y-type distortions by about three 
orders of magnitude. A detection of y ~ 10"* and p ~ 5 x 10"* 
could therefore b ecome possible at the 5a level (see Fig. 12 of 
iKogutetalJboilh , at least in terms of spectral sensitivity. How- 
ever, it is not yet totally clear what the foreground limitations to 
measuring primordial spectral distortions will be. 

Under simplifying assumptions analytic approximation in the 
two extreme cases of p and y-type distortions can be found, how- 
ever, for energy release in the redshift range 10-* < z < few X 10' it 
is clear that a mixture of both types of distortions will form. In that 
case analytic approximations are more difficult and numerical stud- 
ies are required. Furthermore, in the past, mainly the special cases 
of instantaneous energy release at some initia l redshift were numer- 
ically studied in detail for thi s regime (e.g. see lBurigana et al.l 19951 ; 
iBurigana & Salvaterrall2003h . However, from the physical point of 
view cases of continuous energy release over extended periods are 
more natural. In addition, so far only problems with rather large 
energy injection, close to the allowed upper bounds obtained by 
Cobe/Firas, were investigated thoroughly. Again from the obser- 
vational point of view, significantly smaller distortions seem to be 
favoured, and in addition the sensitivities reached by Pixie demand 
studies with very small energy release. To allow forecasting the 
observational possibilities of Pixie with respect to constraints on 
the thermal history of our Universe, significantly refined compu- 
tations of the possible spectral distortions with particular physical 
processes in mind therefore become necessary. 

In this work we revisit the cosmological thermalization prob- 
lem in the light of Pixie. We formulate the thermalization problem 
allowing to resolve very small spectral distortions at levels well 
below the current upper limits of Cobe/Firas. We use improved ap- 
proximations for the double Compton and Bremsstrahlung emis- 
sivities, as well as a detailed treatment o f the cosmological recom - 
bination problem based on CosmoRec jChluba & Thoma3l201 ll) . 
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We explicitly compute the expected CMB spectral distortions for 
different energy injection scenarios. In particular, we discuss two 
scenarios that are even present in the standard cosmological model. 
One is connected with the interaction of CMB photon with adia- 
batically cooling electrons and baryons (see Sect. l2.5.Tt . while the 
other is caused by the dissipation of acoustic waves in the early 
Universe (see Sect. I2.5.2] l. We also demonstrate that the evolution 
of spectral distortions at very low frequencies v < 1 GHz is sig- 
nificantly affected by free-free absorption at late times (z ^ 10^), 
when the Compton process no longer is able to keep the electrons 
and baryons from cooling strongly below the CMB temperature. 
For large energy release this modification is not very significant, 
however, for very small energy release, comparable to ji and y at 
the level ~ 10"**, it becomes important. 

One of the outcomes of this work is CosmoTher^E an open 
source code that solves the cosmological thermalization problem 
for different energy injection scenarios in the limit of small dis- 
tortions. In the future CosmoTherm should be useful for detailed 
forecasts of possibilities to constrain the thermal history of our Uni- 
verse by measuring the CMB energy spectrum. 



2 FORMULATION OF THE THERMALIZATION 
PROBLEM 

A solution of the thermalization problem in the isotropic Universe 
generally can only be obtained numerically. Under simplifying as- 
sumptions several analytic approximations exist, however, these ap- 
proximations have their limitations. For example, for very large en- 
ergy injection, non-linear aspects of the problem become impor- 
tant, while for very small energy release, the precise shape of the 
distortion depends strongly on the small difference in the electron 
and photon temperature, which is comparable to the distortion it- 
self, and is more difficult to account for in analytic approximations. 
Below we provide the required equations for the formulation of the 
problem that allows to include all these aspects consistently. 

2.1 General aspects of the thermalization problem 

To thermalize spectral distortions after some significant energy re- 
lease in the early Universe two main ingredients are needed: (i) the 
number of photons has to be readjusted and (ii) photons have to be 
redistributed over frequency to restore the shape of a blackbody. 

At redshifts z ^ 10^ - lO' electrons and positrons were very 
abundant and electron-positron, electron-electron and positron- 
positron Bremsstrahlung were orders of magnitude more efficient 
in producing soft photons than any other process. Hence thermal- 
ization of CMB spectral distortions was extremely rapid. There- 
fore we shall restrict ourselves to redshifts below z g 5 x lO'. 
By then the number densities of electrons and positrons had frozen 
ouH and only insignificant amounts of positrons were left from the 
era of electron-positron annihilation. During the subsequent evo- 
lution of the (isotropic) Universe the main interactions between 
the photons and baryonic matter are governed by Compton scatter- 
ing (CS), double Compton scattering (DC) and normal electron-ion 
Bremsstrahlung (BR). Furthermore due to the expansion of the Uni- 
verse the photons suffer from redshifting and the non-relativistic 
electrons and baryons cool adiabatically. 



CosmoTherm will be made available at www . Chluba . de/Co smoTherm 
' More precisely the freeze-out of positrons happens at temperatures kT ~ 
m^c-/ 40 orredshift z ~ 5.4 X lO'. 



The photons do not interact directly with the baryonic matter 
(hydrogen, helium and the nuclei of other light elements), but only 
indirectly mediated by electrons, which themselves are strongly 
coupled to the baryons via Coulomb scattering. The timescale, on 
which the electrons and baryons adjust their energy distributions is 
much shorter than any other timescale of importance here. There- 
fore it can be assumed that the electrons and baryons always fol- 
low a (relativistic) Maxwell Boltzmann distribution with one tem- 
perature Te. In this case the problem can be formulated with the 
Boltzmann equations for the evolution of the photon phase space 
distribution function and a coupled equation describing the time 
evolution of the electron (and baryon) temperature. In addition the 
expansion timescale, rexp = 1/^, where H is the Hubble factor, 
is affected by the energy injection process, however, as we argue 
in Sect. |2.4| the corresponding correction is negligible. Below we 
shall write down the photon Boltzmann equation in the expanding 
isotropic Universe and the evolution equation for the electron tem- 
perature with inclusion of heating and cooling. 

2.2 Evolution of the photons in the expanding Universe 

In the early Universe the photons undergo many interactions with 
the electrons in the Universe. As mentioned earlier, the most im- 
portant processes are Compton scattering, double Compton scat- 
tering, Bremsstrahlung and the adiabatic expansion of the Uni- 
verse. Among these processes for most times Compton scattering 
is the fastest. Therefore it is convenient to express all the involved 
timescales in units of tc = l/o"x c, the Thomson scattering time. 

The Boltzmann equation for the evolution of the photon occu- 
pation number, ny, then can be expressed as: 



dny dny 
— H tcv — — 

OT OV 



Any 
"d7 



Any 
d7 



Ajly 

dr 



(1) 



where we introduced the optical depth dr = Atjtc to electron scat- 
tering as the dimensionless time variable. The second term on the 
left hand side is due to the expansion of the Universe and the right 
hand side terms correspond to the physical processes quoted above. 

For computational purposes it is also convenient to trans- 
form to dimensionless frequency. Introducing x = hv/kT-, where 
T- = T-q[1 + z], the expansion term in Eq. (UJ can be absorbed so 
that photons no longer redshift out of the chosen computational 
domain. Here one point is very important: the temperature T.q 
can really be any temperature, as long as it remains constant in 
time. One obv ious choice i s the observed CMB temperature today, 
rn = 2.726 K( lFixsenl2009l) . but we will use this freedom to fix our 
initial conditions conveniently (see Sect. |3.1.Tt . 

For T-,Q = Tq, we have the standard redshift dependent CMB 
temperature, T- = Tqmb = Tq[1 + z]- However, it is important to 
mention that with this definition T, is not necessarily identical to 
the effective (thermodynamic) temperature of the photon field. 



t: 





1/4 


15^3 




= 





1/4 



(2) 



which is found by comparing the total energy density of the photon 
distribution with the one of a pure blackbody. Here ^3 is defined 
by the integral 



&3 



■ J Ax = fff + J x^Any Ax, 



(3) 



over the photon occupation number, . For a blackbody with ther- 
modynamic temperature T- one has ;!pi(x) = l/[e' - 1], and hence 
fff = n'*/l5 a; 6.4939, trivially implying that T* = T-. However, 
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in the presence of spectral distortions after release of energy gen- 
erally r* 9^ T-. This for example means that at higher redshifts 
the effective temperature of the CMB could have been slightly 
lower than today, and only because of energy injection and sub- 
sequent thermalization it has reached its present value. Further- 
more, it is important to mention th at Tp is k nown with precision 
±1 mK iFixsen & Mather 2002: Fi xseijioOSl) . so that deviations of 
the effective temperature from To which are much smaller than this 
are presently indistinguishable. We will return to this aspect of the 
problem again in Sect. l3.1.TI 

In Eq. ([3} we also introduced the distortion, A/i, = ;ji - «pi(x), 
from a blackbody with temperature T-. For numerical computations 
it in general is useful to substitute in this way, as this allows one 
to cancel the dominant terms and thereby linearize the problem, by 
only solving for the (small) correction Awr. 



2.2.1 Compton scattering 

The contribution of Compton scattering to the right hand side of 
the photon Boltz mann equation ([T) c an be treated using the Kom- 
paneets equation jKompaneetsll 195^ : 



d;7v 
dr 



2l a 

dx 



+ (pn^(n^ + 1) 



(4) 



Here we introduced the abbreviations (j> = TJT. and 6. 
This equation describes the redistribution of photons over fre- 
quency, were the effect of Doppler diffusion and boosting, electron 
recoil, and stimulated electron scattering are accounted for. 

The Kompaneets equation was found as lowest order Fokker- 
Planck expansion of the collision term for Compton scat- 
teri ng. We neglec t higher order relativistic corrections (e.g. , 
see [R ephaeli' 199?; Itoh et al. 'l998'; 'Challinor & LasenbiJ 1 19981 ; 
ISazonov & Sunyaey ,1998, ; Nozawa et al, .2006 ). These are ex- 
pected to affect the re sults at the level of a few percent at z > 
10*" - 10' ( IChlubij200f) . so that for the purpose of this paper Eq. igll 
will be sufficient. 

One can easily verify, that the number of photons is conserved 
under Compton scattering by multiplying Eq. (O with x^ and in- 
tegrating over frequency using integration by parts and the fact 
that the photon distribution vanishes sufficiently fast for a" — > 
and X — > oo. Furthermore, one can easily check that for photons, 
which follow a Bose-Einstein distribution, nBE(x) = l/[e'^'' - 1], 
with constant chemical potential fj = po, Eq. ^ vanishes identi- 
cally if (p = 1. For (p i: I in the chosen coordinates the photon 
distribution in kinetic equilibrium with respect to Compton scatter- 
ing is /tBE(x) = l/[e'^^''''^' - 1] with frequency-dependent chemical 
potential fi(x) = fio + x^ - x, where we defined Xe = (px. 

If we now insert n^, = npi(x) + An^ into Eq. ^ and neglect 
terms of (9(An^), we obtain the linear equation 
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1 + e-' 
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(6a) 



(6b) 



1 - 

and A(9e = fe - d-, with e. = k 4.60x 10-'"[1 One can see 
from Eq. ^ that the remaining part arising from the CMB back- 
ground spectrum appears as a source term in the evolution equation, 
which vanishes identically unless T- t T^.. 



Explicitly carrying out part of the derivatives in Eq. ^ and 
reaiTanging terms we find 



dAn,. 



dr 



(9'A;7, 



dx'- 



4 

X 



4 51nf 

X dx 



(7) 



with Z)e = 0e X-, A0 = - 1, and 3, In^ = -2e- V[l - g-^'"]. 

For numerical computations it is convenient to precompute the 
functions (, <f , 5, In f , npj and e"^' once the frequency grid is cho- 
sen. This accelerates the computation significantly. Furthermore, it 
is important to use series expansion of the expressions for x <c 1, 
in particular when differences 1 - e^" are encountered. This is very 
important to achieve numerical stability and the correct limiting be- 
haviour, and will also be crucial for the formulation of the emission 
and absorption term in the next section. 



2.2.2 Double Compton scattering and Bremsstrahlung 

In the cosmological thermalization problem, the double Compton 
process and Bremsstrahlung provide the source and sink of pho- 
tons. Their contribution to photon Bol tzmann equation, Eq. l[l}, can 
be cast into the form ( Rybicki & Lig htmadfl 9791 ; iLightmaiFl 98 ll; 
iDanese & de Zottill982i ; iChlubaii200^ ; 



dra.r 
dr 



= 4[l - n, («"' - 1)] X K(x, d„ 6>e) 



(8) 



where the emission coefficient K is given by the sum of the con- 
tribution due to double Compton scattering and Bremsstrahlung, 
K = A'dc + Kbr- Since K drops off exponentially for x — » cx) 
(see below), the main emission and absorption of photons is oc- 
curring at low frequencies. At small x the photon distribution after 
a very short time is pushed into equilibrium with nl'' = l/[e'' - 1], 
i.e., a blackbody of temperature T- = T^. If once again we insert 
n^ = /ipi(x) -I- A;ji and rearrange terms, we find 



dAii, 



dr 



K{x,9-,9^) 

= 1 ~ ^ f 

DC+BR 1 - e ' x^ 

,A,,a-.-Ox^i-%^, (9) 



with Axe = Xe - X. Note that the first term vanishes if T^. = T-, 
however, the second is pushed to equilibrium only because of inter- 
play with the first term: if a distortion is present then consequently 
Te i= T- and hence the first term either leads to net absorption 
(Je < T-) or emission {T^ > T-) until (at particular frequencies) 
the distortion term is balanced. 



Double Compton scattering 

Due to the large entropy of the Universe at sufficiently high 
redshifts (z > few x 1 0^) DC emission dominates over BR 



jDanese & de Zotti 


19821). The DC emission coefficient 


eiven as jLiehtman 


ll98ll;lThornelll98ll;IChluba et al.ll2007b 



A:Dc(x,a,0e) 



4a 



3n - 



(10) 



where a is the fine structure constant and gdc(x, 6;,6e) is the ef- 
fective DC Gaunt factor. In lowest order of the photon and electron 
energies the DC Gaunt factor factorizes (see_Chluba 2 00S . for more 
details). Furthermore, if the photon distribution is not too far from 
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Figure 1. Effective double Compton conection factor H^^ix). We com- 
pare the result from a full integration of a blackbody spectrum with the 
ap proximation given by E g. il3\ . For comparison also the approximation 
of lBurigana et alj h991bh is shown. Close to the maximum of the CMB 
blackbody spectrum the differences are ~ 20% - 40%. 



full equilibrium one can approximate gdi:(x,d-,dc) using a black- 
body ambient radiati on field and assuming that ~ T,. In this 
case one has (e.g., see lChluball200^ : IChluba et alj2007h 



gi,{x, e., Se) 



1 + 14.16 0; 



(11) 



where if = j x*np\(np] + l)dx = 4n*/l5 ~ 25.976. Here we have 
included the first order relativistic correction in the photon temper- 
ature, however, this term only becomes significant at z > few x 10*. 

The second factor in Eq. jilt allows to go beyond the soft 
photon limit, for which x « 1 was assumed. In lowest order H^^ix) 
only depends on the ambient photon distribution, but is independent 
of the electron temperature. It can be computed using (see Ichlubj 
l2005l for more details) 

HicM « J x'*m,t{x')[l + np,(x' - x)] [^/fo (^)] dx' (12) 

where Ha(w) = [I - 3y + 3y^/2 - y^]/v with y = w[l - w]. The 
factor Hq(w) was first obtained bv loould ( 1 19841) to describe the 
corrections to the DC emissivity when going beyond the soft pho- 
ton limit but assuming resting electrono- In the limit x — > one 

finds wHg(w) — > 1, so that Hi^ix) — » 1. 

Ex pression il2i was also used in the work of lBurigana et al.l 
( Il991bl) . There the approximation H^^ix) x e -""^^^ was given. How- 
ever, as mentioned above with the assumptions leading to Eq. il2i 
the electron temperature is irrelevant, and hence one should set 
0^1. Furthermore, we reexamined the integral and found that 
for background photons that follow a blackbody spectrum 



,3 29 , 11 , 5 4 

1 -I- -X H X- + — .r + — X 

2 24 16 12 



(13) 



provides a much better approximation to the full numerical result 



for Hie (cf. FiglTJ. This approximation was obtained by replac- 
ing npi(x) a e"' and neglecting the induced term in Eq. il2\ . Fur- 
thermore, the resulting expression was rescaled to have the correct 
limit for x — > 0. In particular, for x » 1 Eq. ( 113b captures the cor- 
rect scaling H^^ix) ~ x* e"^" . However, since most of the photons 
are produced at low frequencies x « 1 we do not expect any big 
difference because of this improved approximation. Nevertheless, 
when using the old approximation we found that at early times the 
spectrum is erroneously brought into full equilibrium at very high 
frequencies, just by DC emission and absorption. 

We note here that if the distortions are not small, then in low- 
est order the correction to the DC emission can be accounted for by 
replacing np\ with the solution in the expression for Jjj'. How- 
ever, from the observational point of view it seems unlikely that 
distortions of interest ever exceeded the level A/?v//i,, ~ 10"^, even 
at z ~ 10'. Therefore, the above approximation should be suffi- 
cient. Of course this does not include DC emission from very high 
energy photons that are directly related to the energy injection pro- 
cess. However, in that case the simple approximation used above 
will anyhow need revision, although the total contribution to the 
photon production is still expected to be small. 



Bremsstrahlung 

At lower redshifts (z < few x 10^^) Bremsstrahlung starts to 
become the main source of soft photons . One can define the 
Bremsstrahlung e mission coefficient by (cf. lBuriganaetai]|l991bl : 
lHu&Silkll993al) 



KBR{x,e,) 



X, 0e) . 



(14) 



Note that Hg{w) is 1 /2 of F{w) given by Eq. (27) of lOould h984h . The 
factor of 2 is to avoid double counting of photons. 



2;r Vfei- 4 

Here = h/m^ c is the Compton wavelength of the electron, Zj, 
and gff(Zi,x, 6e) are the charge, the number density and the BR 
Gaunt factor for a nucleus of the atomic speci es i, respectively. Var- 
ious s imple analytical approximations exists l lRvbicki & LightmanI 
Il979h . but nowadays more accurate fitting formulae, valid over 
a wide range of tem perat ures and f requencies, may be found in 
iNozawa etafl l ll998h and lltoh et al] (j2000). In com parison with 
the expressions summarized in Burigana et alj ( Il991bl) we find dif- 
ferences at the level of 10% - 20% for small x. 

In the early Universe only hydrogen and helium contribute 
to the BR Gaunt factor, while the other light elements can be ne- 
glected. In the non-relativistic case the hydrogen and helium Gaunt 
factors are approximately equal, i.e., gn.ff ~ ^He.ff to within a few 
percent. Therefore assuming that the plasma is still fully ionized the 
sum in Eq. l ll4b may be simplified to 2 ~ gn.if N^, where A'b is the 
baryon number density. However, for percent accuracy one should 
take the full expressions for gH.ir and gHe.ff into account, which does 
not lead to any significan t computational burden using the expres- 
sions of lltoh et alj ( I2OOOI) . 

Furthermore, at redshifts z :£ 7000 - 8000 the plasma enters 
the different epochs of recombination. Therefore, the mixture of 
the different species (N^, Hi, Hii, Hei, Heii, Heiii) in the primor- 
dial medium has to be followed. We use the most recent computa- 
tions of the recombination process including previously neglected 
physical corrections to th e recombination dynamics according to 
IChluba & ThomasI ( I2OI ll) . 



2.3 Evolution of the ordinary matter temperature 

As mentioned above, the baryons in the Universe will all follow a 
Maxwell-Boltzmann distribution with a temperature that is equal to 
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the electron temperature. Caused by the Hubble expansion alone, 
the matter (el ectrons plus baryons) temperature would scale as 
Je «: [1 + z]~ dZeldovich et al.ll968h . implying <Ty. However, 
the electrons Compton scatter many times off background photons 
and therefore are pushed very close to the Compton equilibrium 
temperature 



4^3 



T- + 



AJ4 

4^3 



(15) 



within the (distorted) photon field, where J< 



^3 

= / xVv(n.,- + l)dx-If 



the CMB is undistorted, = JjJ /4 and = T., however, in the 
presence of spectral distortions generally T^'^ t T-. Here we also 
defined AJ4 = J4 - Jjj' and Aff^ = ffi - ^3', which for numerical 
purposes is better, as the main terms can be cancelled out. 

In addition, when spectral distortions are present, the matter 
cools/heats via DC and BR, and if some significant energy is re- 
leased, this will in addition heat the medium. Putting all this to- 
gether t he evolution e quation for the electron temperature reads 
(e.g., see IChluball200l for a detailed derivation) 



dpe tcQ 4p. 



4py 



-F = ^ + — -Pe] - — •HDC,BR(Pe) -HtcP,. (16) 
dT QThfy Ch Q'h 

Here we introduced Pe = 1/0 = T^/T-; the heat capacity of the 
mediunO, ka^ = ^klN, + Nh + Nh,] = jkNuU + fn, + X^]- the 
energy injection rate, Q, which for example could be caused by 
some decaying particles (see Sect. 12.51 for more details); and the 
energy density of the photon field in units of electron rest mass, 
Ky0i with Ky = SnAf « 1.760 x 10'"' cm 



Furthermore, 
For numerical purposes 
- 1 , so that this 



Pr = KyOy m 

we defined Pe"* = TTIT- = 1 + ^ - ^ 
it is better to rewrite the above equation in Ape = p, 
formulation for p^^ is very convenient. 

The integral '7Ydc.br arises from the cooling/heating by DC 
and BR, and reads 



r [l -;7v(e* - l)\xK{ 



x,0;,6'e)dx. (17) 



For numerical purposes it again is useful to group terms of similar 
order, when calculating this integral numerically. Usually, 'Hdc.br 
does not contribute very significantly to the total energy balance, 
but for energy conservation over the whole history it is crucial to 
obtain precise values for it. 

2.3.1 Approximate solution for 

Because the timescale for Compton scattering is extremely short 
until z ^ 800, at high redshifts one can solve Eq. \l6i assum- 
ing quasi-stationarity. Even at low redshifts this provides a very 
good first approximation for the correct temperature. Our full nu- 
merical computations confirm this statement. Since the cooling 
from DC and BR is very small, it is furthermore possible to write 

•HDCmip':') « 'WDC,BR(Pe) + d^.'HDC.BRifi'.M^' - p',l SO that 

(18a) 



I + dp.'Hucmip;) + A 
pr * = + - -HdcbrCpI) + •Hdc.br (P^)pi (18b) 



A : 



H fc^h ^ 3 1 + /He + 

4py 8 cr-pXecpy 



H 



(18c) 



' ' We neglected relativistic corrections to th e heat capacity of the electrons, 
which at z,~ lO' would be of order percent jChlubj|200a . 



where p^ is the solution for pe at time-step Tj. Also note that here 
we assumed that the spectral distortion from the time-step t, < t,+i 
is used in all the integrals that have to be taken over the photon 
distribution. This will lead to some small error, which can be con- 
trolled by choosing an appropriate step-size. Furthermore, as we 
will explain below, one can iterate the system at fixed time until the 
solution converges. 

The derivative of "Hdcbr with respect to p^ can be computed 
numerically, however, it is strongly dominated by the term con- 
nected with the exponential factor in the definition Eq. l |17| l. 
For numerical purposes we in general computed all integrals over 
the photon distribution by splitting the spectrum up into the CMB 
blackbody part and the distortion. This allows to achieve very high 
accuracy for the deviations form the equilibrium case, a fact that 
is very important when tiny distortions are being considered. At 
low frequencies we again used appropriate series expansions of the 
frequency dependent functions. 



2.4 Changes in the expansion history and global energetics 

Changes in the thermal history of the Universe generally also im- 
ply changes in the expansion rate. For example, considering a mas- 
sive decaying relic particle, its contribution, px, to the total energy 
density, p,ot, will decrease as time goes by. The released energy is 
then transferred to the medium and because of Compton scattering 
bull{3 of it ends up in the CMB background, potentially leaving a 
spectral distortion, and the matter with an increased temperature. 
Because photons redshift, this implies that the expansion factor 
changes its redshift dependence, and hence the relation between 
proper time and redshift is altered. 

To include this effect, one has to consider the changes in the 
energy density of the photons (or more generally all relativistic 
species) and the relic particles (or more generally the source of the 
released energy). For simplicity we shall assume that the energy is 
released by some non-relativistic massive particle, but this can be 
easily relaxed. We then have 



1 diO'py 
a* dt 
1 daVx 
dt 



-a 



(19a) 
(19b) 



where Py is the energy density of the photons, and a oc [1 + z]"' is 
the scale factor. 

We emphasize that here it was assumed that the matter tem- 
perature re-adjusts, but that the corresponding change in the total 
energy density of the normal baryonic matter is negligible. Further- 
more, Eq. l |19l l implies that all the injected energy ends up in the 
photon distribution. For energy injection at very late times, close to 
the end of recombination, this is certainly not true anymore. In that 
case, it will be important to add the electron temperature equation 
to the system, and then explicitly use the Compton heating term for 
the photons. In this way one has a more accurate description of the 
heat flow between electrons and photons, and only the spectral dis- 
tortion is neglected in the problem. However, in terms of the global 
energetics this should not make any large difference. As we will 
see below, in the case of adiabatically cooling electrons it turns out 



Since the temperature dependent contribution to the total matter energy 
density is negligibly small, well before recombination virtually all the in- 
jected energy will be stored inside the photon distribution. 



© 0000 RAS, MNRAS 000, 000-000 



CMB spectral distortions 1 



that such treatment is necessary in order to define the correct initial 
condition (see Sect. 13. 3t . 

Equation l |19l l has the solutions Py{z) = p^y\z) + ^Py{z) and 



Pxfe) = Pt\z) + Apxfe) with 



Jo 



a\z')Qdz' 
mz')[l+z'] 
Q dz' 



f\z')H(.z')n+z'] 



1 f aHz')Qdz' 



(20a) 
(20b) 



where p^"* = p^.o[l + z]* ~ 0.26 [1 + z]"* eV cm is the unperturbed 
photon energy density, and similarly pi^' = px,o[l + zf- Here pio, 
with i 6 [y, Xj, denotes the present day values of the corresponding 
energy densities. For decaying relic particles with lifetimes shorter 
than the Hubble time this consequently means px.o = 0. We will 
use Eq. l |20at later to compute the appropriate initial condition for 
the energy density of the photon field (see Sect. l3.1.Tt . 

The integrals in Eq. jlOl themselves depend on H(z), however, 
assuming that the total change in the energy density is small, one 
can use H(z) ~ H'-''\z), where //*'''(z) denotes the unperturbed ex- 
pansion factor with 2 = 0. Then the modified expansion factor is 
simply given by H{z) = H''°\z) + A/f(z) with 

AHiz) Apyiz) + Apx(z) 



H(°Kz) 



pT(z) 



- r- 



Q 



dz' 



z pf\z') W<°Hz')[l+z'] 



(21) 



where p|°'(z) = 3 [H'^Xz)]' /SttG is the total energy density of the 
unperturbed Universe. This equation shows clearly that the change 
in the total energy density is arising from the fact that energy re- 
leased by non-relativistic particles is transformed into contributions 
to the photon energy density; if photons had the same adiabatic in- 
dex as matter this would not change anything, a conclusion that 
can be easily reached when thinking about energy conservation. 
However, the additional redshifting of photons implies a small net 
change in the expansion rate. 

In terms of the thermalization problem, the Hubble factor ex- 
plicitly only appears in the evolution equation for the matter tem- 
perature, Eq. l ll6b . while otherwise it is only present implicitly, via 
Tq. For Tq it can be neglected, since this will only lead to an ex- 
tremely small shift in the relation between proper time and red- 
shift that corrects the small spectral distortion, A«i . This conclusion 
can be reached as well when transforming to redshift as time vari- 
able, since then H appears in all terms for the photon Boltzmann 
equation, where AH is a correction-to-correction. A similar reason- 
ing holds for the matter temperature equation. There one could in 
principle replace H with /f"" -I- AH, however, looking at the quasi- 
stationary solution for pe, Eq. l IlSI ). it is evident that this will only 
lead to a minor difference. We checked this statement numerically 
and found the addition of AH to be negligible. We conclude that the 
correction to the Hubble factor enters the problem in second order. 

We also comment that above we made several simplifying as- 
sumptions about the form of the energy release. For example, if 
all the energy is released in form of neutrinos, then both the CMB 
spectrum and the matter temperature remain unaffected, however, 
the Hubble factor will again change, assuming that the energy den- 
sity of the source of the energy release itself does not scale like 
[1 + zT- However, considering this problem in more detail is be- 



yond the scope of this paper and we do not expect major changes 
in our conclusions for this case. 



2.5 Energy injection rates for different processes 

Having a particular physical mechanism for the energy release in 
mind, one can specify the energy injection rate associated with this 
process. Here we assume that most of the energy is going directly 
into heating of the medium, but only very few photons are produced 
around the maximum of the CMB energy spectrum, or reach this 
domain during the evolution of the photon distribution. Also, here 
we do not treat the possible up-scattering of CMB photons by the 
decay products, nor do we allow for Bremsstrahlung emission from 
these particles directly in the regime of the CMB. For a consistent 
treatment of the thermalization problem it will be very interesting 
to consider these aspects in more detail, however, this is far beyond 
the scope of this paper. 

With these simplifications it is sufficient to provide Q in 
Eq. 1 II6I 1 for each of the processes and then solve the thermaliza- 
tion problem assuming full equilibrium initial conditions. Below 
we first discuss the effect of adiabatically cooling electrons and 
the energy injection caused by dissipation of acoustic waves in the 
early Universe. Furthermore, we give Q for annihilating and de- 
caying particles, as well as for short bursts of energy injection. The 
former two processes are present even for the standard cosmologi- 
cal model, while the latter three are subject to strong uncertainties 
and require non-standard extensions of the physical model. 



2.5.1 Adiabatically cooling electrons and baryons 

As mentioned above, if the interactions with the CMB are ne- 
glected, the temperature of the matter scales as Te oc (1 + z)' 
once the baryons and electrons become non-relativistic. At very 
low redshifts {z S 200) it is well known that this difference in the 
adiabatic index of baryonic matter and radiation leads to a large 
differ ence in the CMB and electron temperature jZeldovich et al.l 
Il968l) with Te < Ty, which is also very important for the forma- 
tion ofth££lobal_2j£msignal from the epoch of reionization (e.g., 
see IPritchard & Loebl 120081 and references therein). However, at 
higher redshifts Compton scattering couples the electrons tightly to 
the CMB photons, implying that the baryonic matter in the Uni- 
verse must continuously extract energy from the CMB in order to 
establish ~ Ty until Compton scattering eventually becomes 
very inefficient at z < 200. Consequently, this should lead to a tiny 
spectral distortion of the CMB, tiny because the heat capacity of 
the CMB is extremely large in comparison to the one of matter. 

A simple estimate for this case can be given by equating the 
Compton heating term with the term from the adiabatic Hubble ex- 
pansion: dpyl dt = -Ha\,9- = [m^c^y^ dEj dt. This term is already 
included in Eq. J16I ), so that no extra Q has to be given. With this 
the total extraction of energy from the CMB is roughly given by 



Apy 
Py 



ahOy 
Py.o 



In 



1+z. 



1 +Ze 



a -5.6 X 10"'" In 



1 +Zs 



1 +Ze 



(22) 



For the evolution from z = 2 x 10 until 



10^^ this implies a 



relative energy extraction of 



Icool 



-5.6x10 which is only a 



"cool J 1 1 

few times below the claimed sensitivity of Pixie jKogut et al.l201 ll) . 

Here it is interesting to observe that well before the recombi- 
nation epoch the effective cooling term caused by the electrons has 
the same redshift dependence as annihilating matter (see below), 
i.e. 2 oc (1 + z)*. This implies that the type of distortion will be 
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similar in the two cases, however, the distortion caused by the adia- 
batic cooling of ordinary matter has opposite sign (see Fig.[3j, and 
tend to decrease each other. 

Furthermore, with Eq. l |26t one can estimate the annihilation 
efficiency that is equivalent to the matter cooling term. It turns out 
that if some annihilating relic particle is injecting energy with an 
efficiency /,„„ ~ 1.6 x 10"~''eVs"' then the net distortion intro- 
duced by both processes should practically vanish. As we will show 
below, our computations confirm this estimate (see Fig.[8l(, demon- 
strating the precision of our numerical treatment. 



2.5.2 Dissipation of acoustic waves 

Acoustic waves in the photon-baryo n fluid dissi pate part of their en- 
er gy becaus e of dif fusion dampi ng ( ISilldll968h . As fi rst estimated 
bv lSunvaev & Zeldoyicd lll970ah . and later refined bv lOalvl l ll99lh 
and iBarrow & ColesI ll99lir this effect should lead to a spectral 
distortion in the CMB which depends on the spectral index of 
scalar perturbation, ns, and th erefore could potentially allo w con- 
straining inflationary models l lMukhanov & Chibisovlll"98lh usmg 
future CMB sp ectral measurements jKogut et alJl20T l|). Following 
lHuetalJ l fl994h . the heating rate in a radiation dominated Universe 
can be cast into the form 



(23) 



Ql^ = nn^)PyH(z)[l+z?^"^-'^^\ 

witlOr(x) « 5.2 X 10"' x-^-^ [0.045]'-' for x ~ \. For = 1 this 
implies a total energy dissipation of ~ 3.9 x lO"** between 

z~ 2x10'' and 10^ . This is about 10 times larger than the energy ex- 
tracted in the case of adiabatically cooling matter (Sect. [23Tt over 
the same period. One therefore expects the associated distortion to 
be about an order of magnitude larger. 

Also, like in the case of cooling matter and annihilating par- 
ticles one has Q (\ + z)^ for ;js ~ 1- The equivalent effective 
annihilation rate in this case is ~ 1.5 x 10"^- eV s"'. This sug- 
;ests that the bounds on /ann derived from the CMB anisotro pies 
Slatver et al.ll2009l : lOalli etai]|2009l : iHutsi et alj|200l 1201 ih al- 
ready now might be in tension with this heating mechanism, how- 
ever, the energy release from acoustic waves only goes into heating 
of the medium, while for annihilating particles extra ionizing and 
exciting photons are produced, implying that the effect on the ion- 
ization history could be much smaller here. 

We therefore computed the modification to the recombination 
history caused by this process using CosmoRec, and found a small 
change of AN^/Nc £ 0.06% aro und the maximum of the Thomson 
visibility function jSunvaev & Z eldovich 1970b) close to z ~ 1 100, 
while the electron fraction in the freeze-out tail indeed was affected 
by several percent. The changes in the freeze-out tail of recom- 
bination could have important consequences for the CMB power 
spectra, however, at late times, close to recombination, the simple 
formula Eq. ( 1231 ) should be modified to take into account deviations 
from tight-coupling and also the assumption of radiation domina- 
tion is no longer valid, implying that the rate of energy release is 
overestimated at those epochs. Furthermore, in this case one should 
more carefully treat the fraction of energy that directly goes into 
heating of the matter as opposed to terms appearing directly in the 
photon field. We believe that the overall effect on the recombination 



dynamics arising from this process is negligible. Nevertheless, re- 
fined computations of the energy release rate might become impor- 
tant for precise computations of this problem, but for the purpose 
of this paper the above simplification will be sufficient. 



2.5.3 Annihilating particles 

The rate at which annihilating particles inject energy scales like 
E oc with the number density of the annihilating particles. 
To parametrize the rate of energy injection we use (compare 



lo parametrize the rate or energy iniection ^ 
IPadmanabhan & Finkbei^l2005l : lc"hlubij|2010h 



dE 
d7 



[l-/v]/.nnA'H(l+Z)' 



(24) 



where has dimensions eV s"', and /,, is the fraction of the total 
energy carried away by neutrinos. /„„„ is related to the abundance 
of the annihilating particle, its mass, and the thermally averaged 
annihilation cross section, however, here we treat all these depen- 
dencies with one parameter. 

Current measurements of the CMB anisotro pies already 
tightly constrain r = [1 - /.-]/ann ^ 2 x 10"" eV s"' jsiatver et all 
l2009l : lOalli eTai]|2009l : iHutsi et alj|200^, 1201 ih . This limit is or- 
ders of magnitude stronger than the one obtai ned from Cobe/Firas 
(e.g., see lMcDonald et alJl200lUchluball2010l) . implying that there 
is only little room for CMB spectral distortions from annihilating 
matter, unless a large fraction of the energy is actually released in 
form of neutrinos, i.e. /y ~ 1. 

With Eq l|24j it then follows 



. I ^ gh(z) dE 
m,c- dt 



(25) 



where we parametrized the fr action of the total energy going into 
heating using gh(z). Similar to IPadmanabhan & FinkbeineJ ( |2005|) 

we will adopt guiz) ~ [1 +/He +2(;^p -l-/He(ZHE„ -l-ZHEm)]/3[l -I- /He]. 

Here, respectively Xp, /neZHEu and /hcZheiu denote the free proton, 
singly and doubly ionized helium ions relative to the total number 
of hydrogen nuclei, A'h, in the Universe. At high redshifts, well 
before recombination starts (z > 8000), one has Xj, ~ Zheu ~ 
Zheiu ~ 1, such that gh(z) ~ 1, while during hydrogen recombina- 
tion gh(z) drops rapidly towards gh(z) ~ 1 /3. At low redshifts, a sig- 
nificant part of the released energy is going into excitations and ion- 
izations of ions, rather than pure heating, explaining why gh (z) < 1 
dShull & van Steenberg|ll985l : IChen & Kamionkowskill2004h . 

To estimate the total amount of energy that is injected relative 
to the energy density of the CMB one can simply assume radiation 
domination (l/H ~ 4.79 x 10'''[1 -I- zT' s) and compute the integral 
Apy/py = j ^/pydt. Ill ths CUSS of annihilating matter this yields 



Apy 
Py 



K 5.2 X 10" 



[l-/v],/an 



2 X 10-" eVs-i 



0.75 



0.022 



(26) 



Here we have assumed that /*„„ is independent of redshift. Further- 
more, we assumed that before z ~ 2 x 10'' everything thermalizes 
teurigana et al.|[l991bl) . while effectively no energy is transferred 
to the CMB once recombination ends (z ~ lO"*). 

This estimate showQthat for ~ 2 x lO^^'eV s"' one can 
expect distortions at the level of Aly/Iy ~ 10"'. Interestingly, part of 



We assumed that all cosmologica l parameters asi de from ng have their 
standard values. Furthermore, as in IHu et alj ) 19941) we simply used the Here we neglect corrections because of partial thermalization of distor- 



CoBE normalization for our estimate. tions at ; < 2 x 10 
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the ener gy is released at ve ry high z, where one expects a yU-type dis- 
tortion Ullarionov & Sunyae v 1975a b), however, also at low red- 
shifts some ene rgy is released. This in contra st should lead to a y - 
type distortion jlllarionov & Sunvaev|[l9753lbHHu & Silklll993ah . 
so that in total one expects a mixture of both. We will compute the 
precise shape of the distortion for some cases below (see Sect. |3.5l l, 
confirming these statements. 

2.5.4 Decaying relic particles 

For decaying particles the energy injection rate is only proportional 
to the density of the particle. M otivated by the parametrization of 
IChen & Kamionkowskil J2004l) we therefore write 



time around the heating redshift Zh with some width cth: 



dE 
dT 



: /x Fx A'h e 



-Fx' 



(27) 



Again all the details connected with the particle, its mass and abun- 
dance (here relative to the number of hydrogen atoms in the Uni- 
verse) are parametrized by f^. The exponential factor arises be- 
cause the relic particles rapidly disappear in the decay process once 
the age of the Universe is comparable to its lifetime fx ~ l/fx- 

For particles with very short lifetime a large amount of energy 
could be injected, without violating any of the CMB constraints, 
sin ce at early times the thermalization process is very efficient (e.g. 
see lSunvaev & Zeldovichlll970cl : IHu & SillJll993bl) . This implies 
that in particular early energy release (at z > lO'') could still lead 
to interesting features in the CMB spectrum. Again it will be im- 
portant to see how the characteristics of the distortion change from 
^-type to y-type as the lifetime of the particle increases. 

Under the above circumstances the constraints derived 
from Cobe/Fir as to date are still c omparable to those from 
other probes ijKusakabe et al ] l2006l : iKohri & Takahashil l2010l : 



iKogut et alj|20T l|), at least for particles with intermediate lifetimes 
close to fx ~ 10'- 10'° sec. However, these constraints could be sig- 
nificantly irriprove d with future spectral measurement of the CMB 
jKogut et al.l20T l|). In addition, constraints obtained with the CMB 
temperature and polarization anisotropics strongly limit the possi- 
ble a mount of energy injection by decaying particles wi th long life- 
times dChen & Kamionkowskilliooi : Izhang et al.ll2007l) . 

To estimate the total amount of energy injected by the particle 
we again take the integral Apy/py = / ^ /Py df. This then leads to 



Py 



fi 



8 X 10^ eV 



0.75 



0.022 



1 +zx 



5 X 10-» 



J, (28) 



where for simplicity we once again assumed radiation domination 
(implying t = 1/2H x 2.40 x 10"[1 -I- zY^ s). Furthermore, we 
introduced the redshift corresponding to the lifetime of the particle, 
Zx ~ VTx'o, with to x 2.40 x 10" s. In addition, we defined the 
integral J' = -\f^ e^^ , which for zx ^ is very close to 

unity. Ignoring corrections because of the thermalization procesf^ 
(e.g., seelUu & Silk 1993b; Chluba 2005J, Eq. ^ implies that for 
particles with lifetimes of fx ~ lO's the limits from Cobe/Firas are 
/x ;£ lO^eV for our parametrization. Note that Q\^^^ can be found 
by replacing ^1^^^^,^ in Eq. i25i with the expression Eq. l l27t . 



2.5.5 Quasi-instantaneous energy release 

As next case we discuss the possibility of short bursts of energy re- 
lease. We model this case using a narrow Gaussian in cosmological 



dE 
df" 



^ e-\-'-'h]-a<r\? 



(29) 



V2?!"Crh- 

With this parametrization, assuming sufficiently small (Tj,, one has 



Apy 



« 10" 



2.6 X 10-« eV 



(30) 



We will include these in our estimates of Sect. 13.6.11 



Computationally it is demanding to follow very short bursts of en- 
ergy injection, however, our code allows to treat such cases by set- 
ting the step-size appropriately. Again can be found by replac- 
ing l^i^i^ in Eq. l l25t with the expression Eq. i29\ . 

With the above formula we can also study the effect of go- 
ing from quasi-instantaneous to more extended energy release by 
changing the value of crj,. For short bursts we chose cth ~ 0.05 fh, 
but we also ran cases with cTh ~ 0.25 fh to demonstrate the transition 
to more extended energy release (see Fig.|19t. 



3 RESULTS FOR DIFFERENT THERMAL HISTORIES 

In this section we discuss some numerical results obtained by solv- 
ing the coupled system of equations described above. We first give a 
few details about the new thermalization code, CosmoTherm, which 
we developed for this purpose. It should be possible to skip this sec- 
tion, if one is not interested in computational details. In Sect. l3.3f 
I3.7l we then discuss several physical scenarios and the correspond- 
ing potential CMB spectral distortion. 



3.1 Numerical aspects 

To solve the cosmological thermalization problem we tried two 
approaches. In the first we simplified the computation with re- 
spect to the dependence on the electron temperature. At each time- 
step the solution for the photon distribution is obtained using the 
partial differential equation (PDF) solver devel oped in connection 
with the cosmological recombination problem jChluba & ThomasI 
I2OI ll) . which allows us to setup a non-uniform grid spacing us- 
ing a second order semi-implicit scheme in both time and spatial 
coordinates. However, we assume that the differential equation for 
the electron temperature can be replaced with the quasi-stationary 
approximation, Eq. (TS}, and solved once the spectrum at t,+i is 
obtained. Since both the evolution of the electron temperature and 
the spectral distortions is usually rather slow one can iterate the 
system of equations for fixed r,+i in a predictor-corrector fashion, 
until convergence is reached, and then proceed to the next time. 

In the second a pproa ch, we modified the PDF solver 
of IChluba & ThomasI jioill) to include the additional integro- 
differential equation for the electron temperature. This results in 
a banded matrix for the Jacobian of the system, which has one 
additional off-diagonal row and column because of the integrals 
over the photon field needed for the computation of the electrons 
temperature. Such system can be solved in 0{N log A') operations, 
and hence only leads to a small additional computational burden 
in comparison to the normal PDF problem. We typically chose 
N ~ 4000 - 6000 grid-points logarithmically spaced in frequency 
over the range x ~ 10"^ to 2 x 10^. However, for precise conserva- 
tion of energy we implemented a linear grid at frequencies above 
X = 0.1. Furthermore, we distributed more points in the range 
x = 0. 1 to 15, where most of the total energy in the photon field 
is stored. Even for about 4000 points we achieved very good con- 
servation of energy and photon number. We tested the convergence 
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of the results increasing A' up to 4 - 5 times, as well as widen- 
ing/narrowing the frequency range. 

We then include d the coupled PDE/O DE system into the ODE 
stepper of CosmoRec jChluba et alj2O10l) . This stepper is based on 
an implicit Gear's method up to six order and allows very stable 
solution of the problem in time, with no need to iterate extensively. 
Furthermore, this approach in principle enables us to include the 
full non-linearity of the problem in An^, however, for the formula- 
tion given above this did not make any difference. 

We found that the two approaches described above give ba- 
sically identical results however the second is significantly faster 
and more reliable for short episodes of energy injection. In the cur- 
rent version of CosmoTherm this approach was adopted. For typ- 
ical settings one execution takes a few minutes on a single core, 
however, for convergence tests one execution with A' ~ 20000 and 
small step-size (Az/z ~ lO"*) took several hour. Also for cases with 
quasi-instantaneous energy release a small step-size was necessary 
during the heating phase, such that the runtime was notably longer. 

We also would like to mention that in the current implemen- 
tation of CosmoTherm we assume that the recombination history 
is not affected by the energy injection process and therefore can 
be given by the standard computation carried out with CosmoRec 
This assumption is incorrect at low redshifts for cases with late 
energy release, as small amounts of energy can have a significant 
impact on the ionization history. We plan to consider this aspect of 
the problem in some future work. 



3.1.1 Initial condition 

We performed the computations assuming full equilibrium at the 
starting redshift, Zs- Usually we ran our code starting at z, = 4x 10'. 
However, in order to end up with a blackbody that has an energy 
density corresponding to the measured CMB temperature. To, to- 
day, we had to modify the temperature of the photons and baryons 
at Zs. No matter if the injected energy is fully thermalized or not, the 
energy density of the photon field before the energy injection starts 
should be very close to Py(Zs) ~ Py'"(Zs)-Ap(zs), where Ap(zs) is ap- 
proximately given by Eq. l l20ab . For more precise initial conditions, 
one has to numerically solve the problem for the global energetics 
(see Sect. |2.4t prior to running the thermalization code. 

Although for extremely small energy injections neglecting 
Ap(Zs) t makes a difference in the final effective temperature 
at the endin g redshift, Ze, that is well below the current limit of 
1 mK for To l lFixsen et al.ll20lfl) . for consistent inclusion of the en- 
ergy injection one should compute the initial effective temperature 

of the CMB by r; = T- [l - Ap{z,)/pf\zs)f'* . We then scaled T, 
such that at Zs it is identical to T*. This allows setting A« j = 
which avoids spurious photon production. At the end of the com- 
putation one should therefore expect r*(zc) ~ 7'o[l +Zc] > T^iZs), if 
Q > 0. This furthermore provides a good test for the conservation 
of photons and energy in the code. 



3.1.2 Boundary conditions 

Finally, to close the system, we have to define the boundary condi- 
tions for the photon spectrum at the upper and lower frequencies. 
One reason for us choosing such a wide range over x is that because 
of efficient BR emission and absorption even down to Ze = 200 the 



photon distribution is always in full equilibrium with the electrons 
at the lower boundary. This allows us to explicitly set the spectrum 
to a blackbody with temperature T^. At at the upper boundary we 
also used this condition, which for most of the evolution is fulfilled 
identically, just because of Compton scattering pushing the distri- 
bution into kinetic equilibrium in the very far Wien tail. However, 
at very low redshifts this Dirichlet boundary condition is not fully 
correct, as the timescale on which kinetic equilibrium is reached 
drops. Fortunately, it turns out that this is not leading to any impor- 
tant difference in the solution. 

We confirmed this statement by simply forcing A/;, = at 
both the upper and lower boundary and found that this did not 
alter the solution in the frequency domain of interest to us. Fur- 
thermore, we tried von Neumann boundary conditions of the type 
-I- 0nv(l + "a) = 0. This condition is equivalent to the chosen 
Dirichlet boundary condition, however, the amplitude of the distor- 
tions is left free and only the shape is assumed to be given by a 
blackbody with temperature T^ (i.e., n,- = A(t)/xc for x — > and 
= A(f)e"* for x — > ca). Again this choice for the boundary 
condition did not affect the results well inside the computational 
domain significantly. The current version of the code has both op- 
tions available. 



3.2 Simple analytic description for p- and v-type distortion 

Before looking in more detail at the numerical results, here we 
give a very brief (and even rather crude) analytical description for 
p- and y-type distortions. The formulae pr esented here are moti 
vated by the early works on this problem i Zeldovich & Sunvaev 



'1969: Sunvaev & Zeldovich 1970c; Zeldov ich et al .ll972l ; ISunvaevl 
J974; lUarionov & Sunyaev 1974, 1975a, b), where here we use the 
shapes of the different components as templates to approximate the 
obtained solutions. Some mor e advanced approximations can be 
found in ( iBurigana et alj[l99^ . 

As mentioned in the introduction, for p-type distortions the 
CMB spectrum is described by a Bose-Einstein distribution with 
frequency dependent chemical potential, = l/[e'+''''' - 1]. Pho- 
ton production at very low frequencies always pushed the spectrum 
to a blackbody with temperature of the matter, and at high frequen- 
cies, because of efficient Compton scattering, a constant chemical 
pot ential is found. These limiting cases can be described by ( e . g . , 
see lSunvaev & Zeldovichlll970cl : irilarionov & Sunvaevlll975ah 



P(^) = Po 



(31) 



with the modification that for x « 1 and x » 1 we have the addi- 
tional freedom x + p(x) = x0f, which allows us to renormalize to a 
blackbody with T t T^. Here p„ is the constant chemical potential 
at very large x, and x^ is the critical frequency which depends on 
the efficiency of photon production versus Compton scattering. In 
terms of the brightness temperature this translates into 



T(X): 



(32) 



which both at very large and small x implies T(x) = T-/(f>f. 

The case of y-type spectral distortions is also k nown in con- 
nection with SZ clusters l lZeldovich & SunvaevI 19691) . The spectral 
distortion is given by 



Anl' 



[e- - 1]2 



X -4 

e' - 1 



(33) 



CosmoRec is available at |www . Chluba . de/CosmoRec| 



where is the Compton y-parameter, which directly depends in 
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Figure 2. Evolution of the electron temperature for the standard thermal 
history. The electrons ai'e always slightly cooler than T-. For comparison 
we also show the evolution of the effective temperature of the photon field, 
I - p' = [T. - r*]/r,, which also implies T* < T^. 



the Thomson optical depth and the difference between the elec- 
trons and photon temperature. Again for us only the shape of the 
y-distortion really matters, and we will determine the value of Ve 
from case to case. Most importantly, at small values of x one 
has An;5tV"Pi = ^T/T a -2ye, while at high x it follows that 
A«^°/npi ~ ycX^ and hence AT/T x x. In the case of SZ clusters 
one always has Ve > 0, since the electron temperature is many or- 
ders of magnitude larger than the CMB temperature. However, as 
we will see below (Sect. l3.3T l, in the case of the expanding Universe 
one can also encounter negative values for Ve, even in the standard 
cosmological picture. 

As last component we add a distortion that could be caused by 
the free-free proc ess at very low freque ncies. This can be simply 
approximated by l lZeldovich etaljl972l) 



A"? = >'fr -T 



(34) 



where yg parametrizes the amplitude of the free-free term. 

In the next few sections we present the distortions as a fre- 
quency dependent effective temperature, or brightness temperature, 
which is found by comparing the obtained spectrum at each fre- 
quency with the one of a blackbody: 

T-x 

T{x) = ■ ,. . (35) 

ln(l + n^') 

For approximations of the final distortion we will simply combine 
the three components described above: 



app 

■ 



1 



+ A«;'= + An' 



ff 



(36) 



eA-+p(.v) _ 1 

•f , Ve and yf{ chosen accordingly. If not 



with the parameters x, 
stated otherwise, parameters that are not listed were set to zero. 



3.3 Cooling of CMB photons by electrons in the absence of 
additional energy injection 

As first case we studied the problem of adiabatically cooling elec- 
trons and what kind of distortion this causes in the CMB spectrum. 
As our estimate shows (see Sec. 12.5. ft the distortion is expected to 
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Figure 3. CMB spectral distortion al z = 200 caused by the continuous 
cooling from electrons. Neglecting the spectral distortion in the computa- 
tion of the electron temperature leads to an underestimation of the final 
distortion at low frequencies. We also show a simple analytical fit accord- 
ing to Eq. (36) with parameters /Joo = -2.22 X 10"', = 1.5 X 10"^, 
tpl~ I = -8.0 X Kr'", Ve = -4.3 X 10- "> and Vff = -4.17 X lO^'^. 



be very small. Therefore, this case provides a very good test prob- 
lem for our numerical scheme and we shall see that CosmoTherm 
indeed is capable of dealing with it. 

We started the evolution at = 4 x 10' assuming an initial 
blackbody spectrum (A«,- = 0) with temperature T^-iZs) = T-(zs) = 
r*(zs), where T^Zs) was computed as explained in Sect. l3.1.Tl us- 
ing the approximation Eq. l l20at with Q = -Ha^k Ty for Ap(zs). We 
solved the problem down to Ze = 200, i.e. well after recombination 
ends, however, the distortion froze in before that redshift, so that 
for our purposes Zc = 200 is equivalent to Zc = 0. This statement of 
course ignores other modifications to the CMB spectrum that could 
be introduce d at lower reds hifts, for example by heating because of 
supemovae doh et al. 2003), or shock s caused by large scale struc- 
ture formation HSunvaev & Zeldovicl^l 19721; ICen & Ostriteill 19991; 
iMiniati et alj200(]|) . 



3.3.1 Evolution of the electron temperature 

Figure |2] illustrates the evolution of the electron temperature for 
different settings. The solid line is the result of the full integration 
including all heating and cooling terms, while for the dashed/blue 
line we neglected the cooling of electrons by DC and BR, and set 
Pe'' = 1, meaning that for the temperature equation we ignored the 
spectral distortions that are introduced by this process. As one can 
see, at early times 1 - Pe is about four orders of magnitude larger 
when including all terms than in the case that enforces Pe'' = 1. 
This can be explained by the fact that the heating of the electrons 
actually introduces a small distortion into the CMB, which results 
in a Compton equilibrium temperature Tl°' < T-. At lower redshifts 
the temperatures in the two discussed cases again coincide. This 
is because at those times the cooling by the expansion of the Uni- 
verse starts to dominate over the Compton cooling related to the 
generated spectral distortion. 

For comparison we also show the evolution of the effective 
temperature of the photon field, 1 - p* = [T- - Ty]/T-, which 
also implies T' < T-. We started the computation with T- slightly 
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Figure 4. Evolution of the CMB spectral distortion caused by the continu- 
ous cooling from electrons. At low redshifts one can see the elfect of elec- 
trons starting to cool significantly below the temperature of the photons, 
which leads to strong free-free absorption at very low frequencies. 
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Figure 5. Evolution of the electron temperature for the standard thermal his- 
tory with dissipation of energy by acoustic waves. As a result of the heating 
the electrons are always above T-. For comparison we show the electron 
temperature obtained when neglecting this effect, for which in contrast one 
has Te < T- at all times. 



higher than Tcmb, such that at Ze we correctly have T* ~ Tcmb- 
Numerically, we obtain this value to within ~ 10%, however, when 
computing the initial condition we assumed that the photons were 
cooled all the time. At low redshifts, during the recombination 
epoch this is no longer correct, such that one expects a slightly 
smaller coupling. Furthermore, our computation of the initial con- 
dition assumed no distortion, which again changes the balance to- 
ward slightly lower initial temperature. In particular, the cooling 
electrons start to significantly alter the high frequency tail of the 
CMB distortion (cf. Fig. O, implying a smaller effective temper- 
ature. These aspects of the problem are difficult to included be- 
fore the computation is done. When considering cases in which the 
heating ends well before recombination and is much larger than the 
cooling by adiabatic expansion of the medium, this no longer is a 
problem, since truly bulk of the heat ends up in the photon field. 

We also confirmed this statement by first computing the global 
energy balance problem (see Sect, \2.4i . only neglecting the intro- 
duced distortions. This allowed us to define the initial temperature 
for the run of CosmoTerm more precisely, such that we obtained 
Ty ~ TcMB to within 0.1% at Ze = 200. We conclude that Cos- 
moTherm conserves energy at a level well below 1%. 



3.3.2 Associated spectral distortion 

In Fig. [3] we show the corresponding CMB spectral distortion in the 
two cases discussed above. Here two aspects are very important: 
firstly, the amplitude of the distortion is strongly underestimated 
when one assumes that the Compton equilibrium temperature is just 
T-, i.e., enforces = 1. In this case the distortions do not build 
up in the full way, as the difference of the electron temperature is 
artificially reduced. Since the electron temperature appears in the 
exponential factor of the DC and BR emission and absorption term, 
this leads to a crucial difference at low frequencies. 

Secondly, the distortions at both very low and at very high fre- 
quencies are rather large. This is connected mainly with the low 
redshift evolution of the distortion. Once the Universe enters the 
recombination epochs, the temperature of the electrons can drop 



significantly below the temperature of the photon field (cf. Fig.|2]l. 
This implies significant absorption by BR at low frequencies, and 
also a sizeable down-scattering of CMB photons at high frequen- 
cies, in an attempt to reheat the electrons. Interestingly, the high 
and low frequency distortion is very similar in the two considered 
cases. This also suggests that this part of the distortion is introduced 
at low redshifts, where the electron temperature in both cases is 
practically the same (cf. Fig.|2]l- 

To illustrate this aspect of the problem, in Fig. |4] we present 
a sequence of spectra starting at redshifts during which distortions 
are quickly thermalized (z ~ 10*), passing through the epoch of yu- 
type distortions (z ~ 10^), followed by the y-type era (z ~ 10*), and 
ending well after recombination. Close to the initial time one can 
observe the slightly higher temperature at both low and very high 
frequencies, which is the result of the consistent initial condition. 
At the final redshift the distortion is neither a pure yu- nor a pure 
y-type distortion. At high frequencies it has some characteristics 
of a negative y-type distortion, while around ~ 1 GHz it looks like 
a negative ^-type distortion. At very low frequencies the free-free 
distortion dominates, as explained above. One can see from Fig.|4] 
that the free-free distortion indeed starts to appear at rather late 
times, when the electron temperature departs by more than AT/T ~ 
10"** from the photons. We found that n,. according to Eq. l |36t with 
parameters yu„ = -2.22x 10"', Xc = 1.5xlO--,0f-l = -8.0xlO"'^ 
= -4.3 X 10"'" and yg = -4.17 x 10"'' represents the total 
spectral distortion rather well (cf. Fig.|3j. These effective values for 
fico and ye are several times below the limits that might be achieved 
with Pixie, implying that measuring this effect will be very hard. 

With the values of yu„ and ye one can estimate the amount 
of energy that was release during the /i-era (50000 < z ^ 
2 X 10*) and y-era (z < 50 000), using the simple expressions 
jSunvaev & Zeldovichlll970j) ~ 1.4Apy/py and ye ~ ^Apy/py, 
resulting in Ap-j,/py I,, a 1.6x 10"' and Ap-j,/p-j,|^ « 1.7x 10"'. This is 
consistent with the simple estimates carried out in Sect. 12. 5.71 sup- 
porting the precision of the code regarding energy conservation. 
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Figure 6. CMB spectral distortion at ; = 200 caused by the dissipation of 
energy from acoustic waves for different values of ns . For comparison we 
show the distortion obtained when neglecting this effect. 



Table 1. Spectral distortion caused by the dissipation of acoustic waves. We 
give the parameters for the simple fitting formula, Eq. (36}. For convenience 
we defined Ai^f = (pf - \ . 



ns 


A-oo [10-8] 


xc [10-2] 


A0f [10-8] 


ye[i0-''] 


rff[10-"] 


0.95 


0.735 


1.5 


0.75 


2.2 


1.97 


1.0 


1.75 


1.5 


0.95 


3.5 


2.98 


1.05 


4.05 


1.6 


1.0 


5.9 


4.47 



3.4 Dissipation of energy from acoustic waves 

As next example we computed the distortions arising from the dis- 
sipation of energy in acoustic waves, again starting at = 2 x 10' 
and solving the problem down to Zc = 200. In Fig.|5]we show the 
evolution of the matter temperature and in Fig. [6] we present the 
corresponding spectral distortions in the CMB. In both cases we 
varied the value of the spectral index, ;7s . 

As one can see from Fig. [5] for all cases with dissipation the 
electrons are heated significantly above T-. For < 1 slig htly more 
energ y is dissipated at low z than in the case ns > I (cf . IHu et al.l 
1 1994) . Also the total energy release is smaller in cases with ns < 1- 
This implies that the amplitude of the distortions is larger for larger 
values of iis, and also that the ratio of fi- to v-type contribution 
increases with iis (cf. Table[T). 

As mentioned above, that at late times, during recombination, 
the electrons are significantly hotter when the heating by damp- 
ing acoustic waves is included. We expect that our computation 
significantly overestimates the effect in this epoch, and for a pre- 
cise inclusion of this process during recombination a more detailed 
treatment will be required. In particular the y-type contribution to 
the final spectral distortion could require revision. However, for the 
purpose of this paper the used approximations will suffice. 

Looking at Fig.[6l one can clearly see that in all cases with dis- 
sipation the shape of the distortion is well represented by a yu-type 
distortion around v ~ 1 GHz. At high and intermediate frequencies 
one can again see an admixture of y-type distortion, and at very 
low frequencies, additional free-free emission from electrons with 
r. > Ty becomes visible. We checked that rescaline all the distor- 



Figure 7. Different cosmological signals from the early Universe. 

tions to the same amplitude at y ~ 1 GHz reveals small differences 
in the shape of the distortion in the Pixie bands. However, these 
differences will be more difficult to distinguish, in spite of the total 
amplitude of the distortion being rather large. 

In Table [T] we gave a summary for the parameters needed to 
represent the final distortion in the different cases using the simple 
fitting formula, Eq. l |36l l. The contribution from fi is dominating in 
all cases and is of order yu„ ~ 10-*. By measuring this distortion 
it would in principle be possible to constrain the value of ns in 
addition to the limits that will be obtained with Planck data. Since 
very different systematic effects are involved, this would further 
increase our confidence in the correctness of the inff ationary model. 

Here it is important t o mention that also the cosmologi- 
cal recombination spectru m (lDubrovichlll975l : IChluba & SunvaevI 
l2006l : ISunvaev & Chlubal l2009l) contributes to the signal in the 
Pixie bands. The r esult for the c osmological recombination spec- 
trum computed by IChluba & Sunv aev {200^ for a 100-shell hy- 
drogen atom model is shown in Fig. |7] together with the distortion 
from the damping of acoustic waves. The smaller contributions of 
hehum lOub rovich & Stolvarovlll997l : iRubiiio-JVIartm et alj|20o3: 
IChluba & Sunvae\j|2010h were neglect ed here, however, the pre vi- 
ously omitted free-bound contribution jChluba & Sunvaevl2006l) is 
included. For comparison we also show the distortion arising from 
the interaction of CIVIB photons with adiabatically cooling elec- 
trons alone. Furthermore, we present the total sum of the cosmo- 
logical signals. Here we emphasize again, that even in the standard 
cosmological model all these signals are expected to be present in 
the CIVIB spectrum. 

At high frequencies (v > 700 GHz) the distortion from the 
Lyman-Q' line emitted during hydrogen recombination at z ~ 1400 
dominates. However, in this frequency band the cosmic infrared 
background is many orders of magnitude larger than this signal, so 
that a measurement will be very difficult there. At lower frequen- 
cies, on the other hand, one can see that the cosmological recombi- 
nation radiation leads to a significant contribution to the distortion 
from acoustic damping. It slightly modifies the slope of the dis- 
tortion as well as introduces small frequency-dependent variations. 
These wiggles depend on the dynamics of recombination and there- 
fore in particular should allow to c onstrain the baryon density and 
helium abundance of our Universe tChluba & Sun vaev 2008}. 

From Fig. [7] we can also conclude that when neglecting the 
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effect caused by adiabatically cooling electrons the spectral distor- 
tion arising from the dissipation of acoustic waves is overestimated 
by about ~ 20% at low frequencies. This would lead to a value for 
ng which is biased high. As the amplitude of the distortion is very 
sensitive to the exact value of /is this also demonstrates that refined 
computations of the spectral distortions as well as the heating rates 
because of acoustic damping are required to obtain accurate predic- 
tions in the light of Pixie. Currently, it appears that for ;7s = 0.96 
Pixie could already allow a Icr detection of the effect caused by the 
dissipation of acoustic waves. 

Finally, we would like to compare the result s of our computa - 
tion with the simple analytic predictions given bv lHuetalJh994) . 
There only cases with ns > 1 were discussed, and for iis = 1.05 
they find /j„o ~ 5 x 10"* for a slightly different cosmology, which is 
in very good agreement with the result given here. However, with 
CosmoTherm were also able to compute the admixture of >'-type 
distortions, which at high frequencies dominates. 

3.5 Distortions caused by annihilating particles 

We now focus on the distortions that are created by annihilat- 
ing particles. Recently, in particular the possibility of annihilating 
dark matter with Sommerfe ld enhancement was actively discussed 
in the literature ( see e.g.. iGalH et all |2009|: Isiatver et alj l2009l : 
ICirelli et al.ll2009l ; lHutsi et alj|2009l . l201 iLIZavala et alj|201ll , and 
references therein), resulting in rather tight constraints on this pos- 
sibility from the CMB temperature and polarization anisotropics. 
As mentioned above, these constraints are many orders of mag- 
nitude stronger than those deduced from Cobe/Firas. The reason 
for this is that energy release can directly affect the recombination 
dynamics, depositing energy into the matter, while this amount of 
energy is tiny compared with the huge heat bath of CMB photons. 

Here we therefore take two perspectives: first we compute 
cases that are close to the current upper bound on the annihilation 
efficiency, /ami ~ few x 10""'' eVs"'. Uncertainty in the modelling 
of the effect of dark matter annihilation on cosmological recombi- 
nation can accommodate factors of a few, so that still some range 
is allowed. Second, we compute the distortions for large annihila- 
tion rates to check the consistency of our numerical treatment, and 
illustrate the differences with cases for decaying particles. 

In Fig.[8]we present the spectral distortion after continuous en- 
ergy injection by annihilating particles with different annihilation 
efficiencies, For comparison we also show the result for the 
case without any energy release (see previous section). To under- 
stand aspects of the solution in more detail in Fig.[9]we furthermore 
present the evolution of the electron temperature for two cases with 
annihilation. As Fig. [8] suggests, for ^ 5 x lO^^^eVs"' we 
expect the distortion caused by adiabatically cooling electrons to 
dominate. In this case it will be difficult to improve limits on this 
type of heating process by studying the CMB spectrum. 

For /ann = 10"^' cV s"', ouc Can start seeing a reduction of the 
positive distortion at y ~ I GHz. Also, as discussed in Sect. 13.31 
for /ann = 2 X 10""''eVs"' the net distortion introduced because 
of heating by annihilating particles and the cooling of electrons 
becomes very small distortion. As the annihilation efficiency in- 
creases, the distortion becomes more like a /J-type distortion with 
positive chemical potential at frequencies around y ~ 1 GHz. How- 
ever, at high frequencies we can observe a contribution from a y- 
type distortion with positive >'-parameter. The heating by annihi- 
lating particles at early times pushes the temperature of the elec- 
trons slightly above the temperature of the photons (cf. Fig. O, 
so that these get up-scattered. Depending on the annihilation effi- 
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Figure 8. CMB spectral distortion at ; = 200 after continuous energy injec- 
tion from annihilating particles. 



no annihilation 

4__, = 1 X 10"^' eV s"' 
/ , = 5x 10"eVs"' 




Figure 9. Evolution of the electron temperature for the thermal history with 
annihilating pailicles. Red lines indicate that the electrons are hotter than 
the effective CMB temperature, while blue indicates cooler than this. The 
black/solid line gives the case without annihilation for comparison. 



ciency, this period is more or less extended. For small annihilation 
efficiency this never happens (cf. Fig. |9]l, so that the distortions is 
dominated by negative and );-contributions. 

At low frequencies we can again see the effect of the free-free 
process. Like in the previous section, both the contribution from 
y-type distortion and free-free process are mainly introduced in the 
redshift range corresponding to recombination z ~ 10' - 10'*. In 
particular we also found a characteristic change in the extremely 
low frequency spectrum, when the difference of the electron and 
photon temperature changes sign, however, this occurred at far too 
low frequencies to be worth discussing any further. 

For illustration, in Fig. [TO] we also give some examples with 
very high annihilation efficiencies. As expected the amplitude of 
the distortion increases with /„„■ Furthermore, one can observe a 
shift in the position of the maximal distortion close to y ~ 1 GHz as 
/ann Incrcascs. For annihilation efficiency /^n = 3 x 10"^* eV s"' the 



© 0000 RAS, MNRAS 000, 000-000 



CMB spectral distortions 15 



V [GHz] 




= 5 X 10"^" eV s 
= Sx 10"'' eV s"' 
4__^ = 2x 10"'^ eV s"' 
/ =3x lo'^eVs"' 



Figure 10. CMB spectral distortion aX z = 200 after continuous energy 
injection from annihilating particles with large annihilation cross section. 
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Figure 11. CMB spectral distortion at j = 200 after continuous energy 
injection from annihilating particles in the Pixie bands. 



distortion reaches yu„ ~ 4 x 10"^. This is very similar to the case of 
acoustic damping for «s = 1.05, which would be well within reach 
of Pixie. We also found that the distortion in both cases looks very 
similar, however, current limits on /an„ from the CMB anisotropies 
disfavour so high annihilation efficiencies. 

Finally, in Fig.[TT]we show the distortions for some of the pre- 
vious cases, but focused on the spectral bands of Pixie. One can 
see that in all shown cases the distortions to lowest order have a 
very similar shape, while only the amplitude is varying from case 
to case. Furthermore, the typical amplitude of the distortions is 
hT IT ~ 10"' - 10"^, rendering a measurement in terms of the 
sensitivity of Pixie rather difficult. At frequencies below 30 GHz 
the situation would be slightly better, as the distortion increases 
to the level of AT/T > 10"**. In this fre quency band also the dis- 
tortions from the recomb ination epoch jChluba & Sunvaevll2006l : 
ISunvaev & Chluball2009h would become visible (see Fig.|7). 

Presently, we conclude that improving the limits on energy 
injection by annihilating particles using Pixie will be very hard 



unless /a„n exceeds ~ 10^'^eVs"'. With the current upper limits 
on the effective annihilation rate is also seems that in comparison 
to the distortion introduced by the dissipation of acoustic waves, 
those arising from annihilation would only contribute at the level 
of 10% - 20%. If annihilating particles are present in the Universe 
this might lead to problems in interpreting future spectral measure- 
ments that could be carried out by Pixie. This is especially severe, 
since the shape of the distortion in both cases is very similar, im- 
plying that more detailed forecasts will be required to assess the 
observational possibilities with Pixie, also including possible ob- 
stacles introduced by foregrounds (see comments below). 

3.6 Distortions caused by decaying particles 

In the case of decaying particles, observational constraints are less 
tight. Therefore one can in principle allow large amounts of energy 
injection well before the recombination epoch. However, here we 
are not so much concerned with deriving detailed constraints on 
this process using our computations. We rather wish to demonstrate 
how the shape of the distortions depends in the particle lifetimes. 
To this end we therefore chose parameter s which are not in tensio n 
with the limits deduced from Cobe/Firas jpixsen & Mathedl2002l) . 

In Fig.[T2]we compiled the results obtained for fixed total en- 
ergy release by particles of different lifetimes. One can clearly see 
that for particles with short lifetimes, corresponding to zx ^ 2x10*, 
only a small part of the released energy remains visible as spec- 
tral distortion. This is because the processes of thermalization is 
very efficient, and does not allow much o f the distortion to survive 
teurigana et al.ll99Ibl : lH"u & Silkll993l. cf. also). As zx decreases 
to Zx ~ 5 X lO"* the low frequency feature at v ~ 1 GHz becomes 
more visible. The overall shape of the distortion is well character- 
ized by a yu-type distortion, and only at very low frequencies it de- 
viates slightly, because of the free-free distortion introduced at late 
times. Since in the chosen example the effective energy release rate 
of the decaying particles is rather large, the free-free distortion is 
not as visible as in the cases discussed in the previous two sections. 

As we decrease zx down to zx ~ 5 x 10"* we can see a change in 
the characteristic spectral distortion. The low frequency /J-type fea- 
ture starts to disappear, while at intermediate and high frequencies 
a significant y-type distortion starts to mix in. Decreasing zx even 
further the distortion clearly starts to be dominated by a y-type dis- 
tortion, with a flat AT/T at intermediate frequencies and a charac- 
teristic rise of AT/T at high frequencies. However, at low frequen- 
cies the interplay between y-type and free-free distortion becomes 
important, leading to another positive feature at v ~ 500 MHz. 

To understand the effect of decaying particles on the CMB 
spectrum a little better in Fig. [T3] we present the evolution of the 
electron temperature for some cases of Fig. [12] One can see that 
for decreasing values of zx at high redshifts the electron tempera- 
ture follows the case without energy injection for a longer period. 
Then, once the heating by decaying particles becomes significant, 
the electron temperature becomes larger than T-. After the heating 
stops for cases with zx ^ 10^ the relative difference in the elec- 
tron temperature remains rather constant, with only slow evolution. 
Caused by the heating the effective temperature of the CMB also 
increased and after it ceased the electrons simply keep the temper- 
ature dictated by the distorted CMB photon field. 

In cases with zx :£ 10^, however, one can observe an extended 
period after the maximal heating at which the electrons loose some 
of their heat again. Having a closer look at the cases with zx ^ 10' 
one can find the same there, but much less pronounced. At high 
redshifts the Compton interaction is extremely fast and allows the 



© 0000 RAS, MNRAS 000, 000-000 



16 Chluba and Sunyaev 



V [GHz] 
10 




V [GHz] 
10 



= 5x10" 
= 2.5 X lo' 
= 1 X lo' 
= 7.5 X 10^ 



V [GHz] 
10 




Figure 12. CMB spectral distortion at z = 200 after energy injection from 
decaying relic particles. In all cases we fixed /jec = 2zxeV, which cor- 
responds to a total energy release of Apy/p 



1.3 X 10"". For the 

effective temperature of the CMB this implies AT*/Tcmb 3.2 X 10"^ at 

Zs = 2x lO'' and at Ze = 200 in all cases we found \AT'/Tcmb\ ~ 10"'". 




Figure 13. Evolution of the electron temperature for the thermal history 
with decaying particles. Parameters were chosen like in Fig. 1121 Red lines 
indicate that the electrons are hotter than the effective CMB temperature, 
while blue indicates cooler than this. The black/solid line gives the case 
without annihilation for comparison. 



temperatures of electrons and photons to depart only slight, even 
with significant energy release. At low redshifts, Compton scat- 
tering becomes much less efficient, so that during energy release 
larger differences between electrons and photons are possible. Dur- 
ing these periods the electrons are notably hotter than the CMB, so 
that photons get up-scattered and a >'-type signature can arise. 

In Fig. [14] we illustrate the evolution of the CMB spectral dis- 
tortion caused by the heating from decaying particles with diff'er- 
ent lifetimes. The upper panel gives an example for a particle with 
short lifetime. The distortion clearly is close to a //-type distortion 
until very late times. The only difi'erence is because of the effect 
of electrons cooling significantly below the CMB temperature at 
late times, introducing a small modification because of free-free 
absorption in the 100 MHz frequency band. In the central panel we 
give a case which at the end has the character of both yu and y- 
type distortion. Initially, it starts like a yu-type distortion, but heating 
continues to be significant down to z ~ 10', when electrons obtain 
temperatures larger than the CMB, such that photons are partially 
up-scattered. At the end of the evolution the spectrum remains in a 
state that is a mixture. Finally, in the lower panel of Fig.[T4]we give 
an example for a case that looks like a pure y-distortion at high and 
intermediate frequencies. In this case, energy is mainly released at 
times when Compton scattering is unable to re-establish full ki- 
netic equilibrium with the electrons. However, at low frequencies 
one can again observe the effect of cooling electrons during at after 
the epoch of recombination. 

In Fig. [TS] we show the distortions for some of the previous 
cases, but focused on the spectral bands of Pixie. In contrast to the 
case of annihilating particles, where the shape of the distortion was 
rather insensitive to the effective annihilation rate, for decaying par- 
ticles the shape of the distortion varies strongly with its lifetime. 
This should make it possible to distinguish the effect of decaying 
particles from the other sources of energy release discussed so far. 
For the chosen energy injection rate the typical amplitude of the 
distortions is AT/T ~ 10"' - 10"^, which is well within reach of 
the Pixie sensitivities. However, to forecast the possible constraints 
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Figure 14. Evolution of the CMB spectral distortion caused by the heating 
from decaying particles with different lifetimes. At ; < lO** one can see the 
effect of electrons starting to cool significantly below the temperature of the 
photons, which leads to free-free absorption at very low frequencies. 




Figure 15. CMB spectral distortion aX z = 200 after energy injection from 
decaying relic particles with different lifetimes in the Pixie bands. In all 
cases we fixed /dec = 2 zx eV, which con'esponds to a total energy release 
of Apy/p. 



1.3 X 10"°. For the effective temperature of the CMB this 

implies Ar*/TcMB 3.2 X 10"'' at = 2 X lO' and at Ze = 200 in all 

cases we found IAT'/TcmbI ~ 10"'". 



from Pixie requires consideration of more cases and realistic fore- 
ground models. 



3.6.1 Upper limits from analytic estimates 

lHu&Silkiri993bh provided simple analytic expressions that allow 
estimating the final spectral distortion after some energy release 
caused by decaying relic particles. These expressions were widely 
used in the literature to place limits on the possible amount of de- 
caying particles with different lifetimes, and here we wish to com- 
pare them with the results of our computations. 

To obtain the analytic estimates one can start with the simple 
approximations for single energy release at Zh- At high redshifts, 
during the era of ^-type distortions o ne haj"l l lDanese & de ZottH 
ll982l : lHu&s"ii3ll993al : IChluball2005[) 



where the thermalization redshift 
Zu = 1.98 X lO*- 



-([l+Ch]/[l+:;,])^'^ 



"i-yp/2 


-2/5 




-2/5 


0.88 




0.022 





(37) 



(38) 



was already used several times above. For tih ISunvaev & ZeldovichI 
l ll970cl) gave the well known approximation p^ « 1.4 — , where it 
is assumed that a negligible amount of photons is injected, but bulk 
of the energy comes out as heat. 

In Eq. l |39l l, the exponential factor acts as a visibility function 
for spectral distortions. At redshifts z < practically all energy 
ends up as CMB spectral distortion, while at z > thermalization 
exponentially suppresses the residual distortion with double Comp- 
ton emission being the main source of photons. 



" In a baryon dominated Universe BR is more important than DC emis- 
sion. In this case the one has a /i,, e"(['+=hl/[l+:brl)''"' withzbr ~ 6.2x10* 
iSunvaev & Zeldovichlll970cl) . 
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Figure 16. Resulting efficiency integral J' for decaying particles with life- 
times corresponding to redshift zx- 



To compute the total distortion arising in the yu-era from de- 
caying particles one simply has to calculate the spectral visibility 
weighted energy release rate: 



Py 



^ r 1 d£ 



„-([l+c]/[l+c^]V'^'- 



(39) 



which, assuming radiation domination, in our parametrization, 
Eq. ( 127 1 reads 



Apy 

Py 



fi 



8x lO^eV 
where we defined the integral 



0.75 



0.022 



1 +zx 



5x 10* 



Jo 



5/2^-5/4] 



-0r 

- 211/18 j4/9 ^10/9 



exp 



9A 



10/9 



28/9 55/9 



J, (40) 

(41a) 
(41b) 



with Ax = [I + Zx]/[1 + Zf,]- The integral J can be easily solves 
numerically and is shown in Fig. [16] together with the result from 
the approximation, Eq. l|41bt . which works very well for zx ^ Zf,- 

In the work o f lHu & Silk (1993hl) , this estimate was performed 
in a slightly different way. There it was assumed that all the en- 
ergy released by the decaying particles effectively is injected at 
time feff ~ fx. This can be concluded from Eq. (8) in their pa- 
per, where the exponential factor reads e"*'iic/'x)''* = g-'^x , which 
implies Jhu a e x . In Fig. [16] we also plotted this version for 
jf and find that for zx > z^ it stron gly unde r estima tes the actual 
value of J^, as already pointed out bv lChlub^ ( l2005l) . This implies 
that the limits derived from Cobe/Firas for particles with lifetimes 
tx S 6 X 10* s are significantly stronger. 

Numerically we were able to compute the efficiency function 
jf using CosmoTherm. In practice just defines how much of the 
energy that was released remains visible as spectral distortion to- 
day. Assuming a constant total energy release, one can therefore 
compute simply varying the lifetime of the particle and com- 
paring the effective value for /j„ with the total amount of injected 
energy. To make the results more comparable we switched off BR, 
since for the estimate above this was not included consistently. This 











/ 4 = 5xl0' 


late free-free 1- / .-y 


- 


~ absorption \; ■! 


- 




no Bremsstrahlung 




DC&BR 




analytic approximation 



Figure 17. CMB spectral distortion at z = 200 after energy release by de- 
caying particles with redshift zx = /x = 5 X 10* and Ap/p ~ 6.3 X 10"^. 
In this figure the importance of BR is illustrated. Also we show simple an- 
alytic approximations according to Eq. 1321 with poo = 2.65 X 10"** and 
x= 1.6 X 10"^ in the case with BR, andp„ = 1.8 X 10"* andx = 1.6 X 10"^ 
in the other. 



also makes it easier to define p„ as the late changes in the distortion 
at low frequencies are not arising (see Fig.|17|l. From Fig.[T7]we can 
also see how much the low frequency spectrum is affected by BR. 
The position of the maximal temperature dip in the case without 
BR is close to x ~ 1.6 x 10 ^ while with BR it is at x ~ 1.6 x 10"^. 
This demonstrates t he well known fact that D C becomes inefficient 
at low redshifts (see lPanese & de Zotti]ll982h . 

The result of this exercise is also shown in Fig. [16] for 
Apy/py ~ 6.4 X 10"*. As one can see, the agreement with the ana- 
lytic estimate is excellent for this amount of energy injection. How- 
ever, for larger energy injection we found that < at z > z^. 
Also, when switching on BR the simple analytic formula cannot be 
directly applied, as in those cases the distortion no longer is just a 
simple p-type distortion, but the low redshift and frequencies other 
contributions can be noticed. In these cases, full numerical compu- 
tations for each case should be carried out. 

We also comment here that, using the simple expression, 
Eq. l l41bt . it is clear that the distortions introduced during the 
epoch of electron-positron annihilation (z ~ 6 x 10') are com- 
pletely unobservable. Although the total energy release from this 
epoch is Ap/p ~ 1, the distortions are suppressed by at least 
J ~ 5 X 10"^°. For this simple estimate the very efficient electron- 
electron and electron-positron Bremsstrahlung process have not 
been included but would lead to even faster thermalization. Fur- 
thermore, e ven only with normal BR on e would reach the same 
conclusion dSunvaev & Zeldovichl[T97o3) . Anything happening at 
z ^ 10' with Apy/py < 10"^ will lead to distortions that are no 
larger than p„ ~ 10"'. 



3.7 Distortions caused by quasi-instantaneous energy release 

To close our computations of spectral distortion after energy re- 
lease we considered some examples with quasi-instantaneous en- 
ergy injection. The evolution of the electron temperature is reported 
in Fig. [T8] and a few cases for the final spectral distortions in the 
Pixie bands are shown in Fig. [19] We started your computation 
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Figure 18. Evolution of the electron temperature for thermal histories with 
quasi-instantaneous energy release. The total energy release was Ap/p ~ 
10"^ in all cases. The black/solid line gives the case without energy injec- 
tion for comparison. Red colour indicates that T,. > T- while blue colour 
means <T-. 




x 



Figure 19. CMB spectral distortion at z = 200 after quasi-instantaneous 
energy injection at different redshifts in the Pixie bands. The total energy 
release was Ap/p ~ W-^ in all cases. The first four curves were computed 
for o"!, = 0.05 while for the last one we used tri, = 0.25 ft,. 

at z = 4 X 10^ with the initial condition defined as explained in 
Sect. [3TT] 

The electron temperature rises steeply close to Zh and falls off 
very fast once the energy release ceases. For the spectral distortion 
one can again observe the diflterence in the characteristic spectrum, 
which is dominated by /J-type contributions for cases with large, 
and is clearly of y- type for very small values of Zh- For zi, ~ 5 x lO"* 
the distortion is a mixture of both. 

The differences in the shape of the distortions are very visible, 
however, comparing to the case with decaying particles shows that 
quasi-instantaneous energy release gives rise to very similar distor- 
tions for injection at equivalent epochs. This suggests that distin- 
guishing these two cases could be rather demanding. 



To make this point even clearer, for Zh = 5 x lO"* we also 
ran a case with cr\, = 0.25 /h and the same total energy injection. 
The result is also shown in Fig. [19] Clearly, it is rather hard to 
distinguish the two lines corresponding to Zh = 5 x 10', however, 
the differences are still at the level of a few percent. For sufficiently 
large energy release this could possibly be measured with Pixie, 
however, more detailed forecasts will be necessary. 



4 DISCUSSION AND CONCLUSION 

We presented detailed computations of the CMB spectral distor- 
tions introduced by different physical processes. We focused in par- 
ticular on small distortions that could be within reach of Pixie. For 
this purpose CosmoTherm, a thermalization code which solves the 
couple Boltzmann equation for photons and electrons in the ex- 
panding isotropic Universe, was developed. Improved approxima- 
tions for the double Compton and Bremsstrahlung emissivities, as 
well as the latest treatment of the cosmological recombination pro- 
cess were taken into account. 

We demonstrated that the interaction of CMB photons with 
adiabatically cooling electrons and baryons results in a mixture of 

a negative fi- and y-type distortion with effective yu 2.2 x 10"' 

and y 4.3 x 10~'° (see Fig.[3] l. For the currently quoted sensi- 
tivities of Pixie l lKogut et al.l201 ih this effect is unobservable, even 
in the most optimistic case of no foregrounds. However, the rapid 
progress in detector technology and the possibility of extending the 
frequency bands of Pixie below 30 GHz, might render this process 
interesting in the future. 

Furthermore, we computed the shape of the distortion arising 
from the dissipation of acoustic waves in the expanding Universe 
(se e Fig. [Sjl using estimates for the energy release rates provided 
bv lHuetalJh994h . We find an effective yu ~ 8.0 x lO"' with an 
admixture of y ~ 2.5 x 10"' for «s = 0.96. It therefore appears that 
for Hs = 0.96 Pixie could already allow a lo" detection of the effect 
caused by the dissipation of acoustic waves. 

Nevertheless, the distortion created by the damping of acous- 
tic waves might be precisely measured in the future. This would 
in principle allow us to place additional constraints on different in- 
flationary models, however, as we discuss here for precise predic- 
tions of the expected distortion one should also include the contri- 
butio n from the cosmological recombination rad iation ( Du broviciil 
1 19751 : IChluba & SunvaevI |2006|; ISunvaev & Chluba 2009.) . which 
was emitted at z ~ 1400, and contributes significantly at very high 
and low frequencies (see Fig.|7j. In addition, our computations in- 
dicate that the energy injection rate associated with the dissipative 
heating process during the recombination epoch is very likely over- 
estimated and should be refined for reliable forecasts. This could in 
particular affect the admixture of y-type distortions arising from 
this heating mechanism, and we plan to look at this problem in a 
future publication. We also emphasize that all the aforementioned 
distortions should arise in the standard cosmological model. Ob- 
serving the cosmological recombination spectrum together with 
these distortions might therefore open a way to place additional 
constraints on the Universe we live in. 

We also computed the possible distortions arising from an- 
nihilating and decaying relic particles. In the case of annihilating 
particles the distortion is dominated by a ^-type contribution (see 
Fig. [U, however, at very low frequencies the free-free processes 
leads to significant modifications. Similarly, at high frequencies a 
late y-type contribution arises. The shape of the distortion is very 
similar for different annihilation efficiencies, but the amplitude and 
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sign depend significantly on this. For very small annihilation rates 
the distortion arising from cooling electrons dominate and make 
the annihilation signal very hard to observe. Furthermore, if an- 
nihilating particles are present in the Universe, their effect on the 
CMB spectrum could be confused with the one from dissipation of 
acoustic waves, however, for current upper limits on the effective 
annihilation rate, this is likely restricted to the level of 10% - 20%. 
We furthermore conclude that presently improving the limits on en- 
ergy injection by annihilating particles using Pixie will be very hard 
unless /a„n exceed ~ 10"-^ eV s"' . 

In the case of decaying particles, the shape of the distor- 
tion strongly depends on the lifetime of the particle (see Fig. 1151). 
This should allow distinguishing between different particle mod- 
els rather easily, unless the energy deposition rates become too 
small. We also revisited the simple estimates for the amplitude of 
the spectral distortion in the yu-era, showing that the residual dis- 
tortion caused by decaying particles with lifetimes shorter than 
fx ~ 6 X 10* sec was underestimated. This implies that the CMB 
limits on decaying particles are possibly slightly stronger than pre- 
viously anticipated (see Fig.ll6t. However, since the distortions for 
particles with lifetimes fx ~ 10* - 10' sec are neither pure /j- nor 
pure y-type deduced upper limits are more difficult to interpret. 

Finally, we also computed the distortions arising in cases with 
quasi-instantaneous energy release (see Fig.[T9}- These distortions 
are very similar to those obtained in cases with decaying parti- 
cles, so that distinguishing the two mechanisms will be challeng- 
ing. However, a detailed forecast of the observational possibilities 
with Pixie regarding all these different cases should be performed 
to reach a final conclusion about this. These forecasts should in- 
clude the possibilities of simultaneous energy release mechanisms. 
In addition, the distortions introduced a t lower redshift s, for exam- 
ple by heating because of supernovae jOh et al.ll2003l), or shocks 
during large scale structure formation jSunvaev & Zeldovichll97^ : 
ICen&Ostrik"e3ll999l :l Miniati et al. I l2000h . have to be considered. 
Also, the effect of unresolved SZ clusters jMarkevitch et al.lll991 ). 
the k inetic SZ effect from reionization (e.g., see iMcOuinn et al.l 
and references therein), and the integrated signals from d usty 
galax i es (e.g., se e Righi et al. 2008; Viero et al. 2009; Vieira et al.l 
I2OIOI ; iDunklev et alj|201(]| ; IPlanck Collaboration et alj|201 ih will 
contribute at an important level. In all these cases one should con- 
sider the possibilities to include spacial-spectral and polarization 
information to separate the different components. However, we 
leave a detailed investigation for some future work. 
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